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Static spherically symmetric perfect fluid solutions are studied in metric f(R) theories of gravity. 
We show that pressure and density do not uniquely determine f(R) ie. given a matter distribution 
and an equation state, one cannot determine the functional form of f(R). However, we also show 
that matching the outside Schwarzschild-de Sitter-metric to the metric inside the mass distribution 
leads to additional constraints that severely limit the allowed fluid configurations. 
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I. INTRODUCTION 

Modern day cosmological observations such as super- 
novae type la cosmic microwave background [H and 
large scale structure [3| provide strong evidence against 
a critical density matter dominated universe. Instead, 
the current cosmological concordance model is a criti- 
cal density universe dominated by cold dark matter and 
dark energy in the form of some kind effective cosmolog- 
ical constant. The traditional cosmological constant is 
maybe the leading dark energy candidate (for a review 
see e.g. 0]), but a large number of other alternatives 
have been studied in the vast literature on dark energy. 

Modifying general relativity (GR) to explain the 
present day acceleration is an often considered avenue 
of research. In particular, a class of models that has 
been extensively studied in the recent years are the f(R) 
gravity models that replace the Einstein-Hilbcrt (En- 
action with an arbitrary function of the curvature scalar 
(see e.g. 0, i, 0, 0, U M, HI E Q and references 
therein) . Modifying the gravitational action is faced with 
many challenges, however, and obstacles such as insta- 
bilities [HI, EE EH as well as constraints arising from 
known pro perties of gravity in our solar system (see e.g. 
[H, l20Ll2l| and references therein) need to be overcome. 
Also the large scale perturbations present a challenge for 
f(R) gravity theories [H, HI]. One of the most direct 
and strictest constraints on any modification of gravity 
comes from observations of our solar system. This is of- 
ten done conformally transforming the theory to a scalar- 
tensor theory and then conside ring the Parameterized 
Post-Newtonian (PPN) limit [H,[2J (see also [H,[26j] for 
a discussion). The question of Solar System constraints 
on /(i?)-theories has recently been extensively discussed 
by a number of authors [29L T30L l3lL l3^| . 

The essence of the discussion is the validity of the 
Schwarzschild-de Sitter (SdS) -solution in the Solar Sys- 
tem. The SdS-metric is a vacuum solution to a large class 
of /(-R)-theories of gravity. However, due to the higher- 
derivative nature of the metric /(-R)-theories, it is not 
unique and other solutions can also be constructed (see 



eg. [15, 39]). This fact is also present in the cosmological 
setting, rendering any cosmological solution non-unique 
and hence the form of f(R) cannot be uniquely deter- 
mined from the expansion history of the universe alone 

In the recent literature this question has now been ad- 
dressed without resorting to scalar- tensor theory [29|, [3(| 
[3lj] . The result is compatible with the scalar-tensor the- 
ory calculations: the Solar System constraints are valid 
in a particular limit that corresponds to the limit of light 
effective scalar in the equivalent scalar-tensor theory. In 
terms of the /(i?)-theory, this is equivalent to requiring 
that one can approximate the trace of the field equations 
by Laplace's equation [3(|. As a result, the often consid- 
ered 1 /R theory is not consistent with the Solar System 
constraints in this limit, if the 1/R term is to drive late 
time cosmological acceleration. 

In this work, we approach the question differently by 
asking what kind of a mass distribution is required so that 
the SdS-metric is the solution outside a spherically sym- 
metric body. In other words, given the SdS-metric, what 
is the matter distribution that has the correct boundary 
conditions. Our approach is general, we do not make any 
assumptions about the /(i?)-theory. 

Here we show that like in the cosmological setting, the 
mass distribution alone cannot in general determine the 
gravitational theory, or the functional form of f{R). Im- 
posing the SdS-metric as a boundary condition docs limit 
the allowed solutions however. We also give a prescrip- 
tion how one can in principle solve the mass distribution 
that has the SdS-metric as the outside solution, given a 
gravitational theory and the equation of state of matter. 



II. f(R) GRAVITY FORMALISM 

The action for f(R) gravity is (8irG = 1) 



S = 



r ~g[f{R) + C, 



(1) 



The field equations resulting in the so-called metric ap- 
proach are reached by variating with respect to g^: 
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F(R)R pi/ - -f(R)g^ - V t y i/ F{R)+g^UF{R) = T™, 

(2) 
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where T™ is the standard minimally coupled stress- 
energy tensor and F(R) = df/dR. Alternatively, this 
can be written in a form similar to the field equations of 
General Relativity (GR): 



1 



(3) 



where the stress-energy tensor of the gravitational fluid 
is 



T, 



l r l 



+ F(Ry a P(g atl g Pu ~g^g a ^} (4) 



and we have defined 



T™=T™/F{R). 



(5) 



Contracting Eq. ([4| and assuming that we can describe 
the stress-energy tensor with a perfect fluid, we get 

F{R)R - 2f{R) + SDF(R) = p - 3p. (6) 



III. SPHERICALLY SYMMETRIC PERFECT 
FLUID SOLUTIONS 

We are interested in static, spherically symmetric per- 
fect fluid solutions, or as commonly referred to in the 
literature, SSSPF-solutions. Study of such solutions can 
be traced back to Schwarzschild [36[ and the literature 
is extensive (see eg. (37j for reviews). The surface of 
the fluid sphere is defined by the surface of zero pressure 
where the interior solution is matched to the outside met- 
ric. 

In spherically symmetric coordinates the metric can 
generally be written as 
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(7) 



The corresponding continuity equation is 



P'(r) lA'jr) 
p(r)+p(r) 2A(r)' 

where prime indicates a derivation with re 
d/dr. 



respect to r, ' = 



A. Uniqueness 

The higher derivative nature of metric f(R) theories 
of gravity can lead to non-uniqueness of solutions of the 
field equations. For example, the cosmological solution 
a(t) for a given f(R) theory is not unique but other f(R) 
theories with the same expansion history exists [22 ] . This 



also true for vacuum solutions: even though the SdS met- 
ric is a solution of the field equations in vacuum, others 
also exist [39j . 

Here we consider the question of uniqueness in the 
presence of matter. This is most conveniently done by us- 
ing a form of the field equations where f(R) is eliminated 
from the equations by using the contracted equation Eq. 



FRuv - ~ (FR-UF)g^ - V M V„F - T™ - I {p-3p)g^. 

(9) 

Given a matter distribution, p(r), p(r), one can solve 
for A from the continuity equation and substitute into 
Eq. © giving a set of differential equations for F(r) and 
B{r). 

This set can be solved algebraically for B(r) and B'(r) 
so that both B(r) and B'(r) can be expressed in terms 
of F(r) 7 p(r), p(r) and their derivatives (see Appendix 
for details). Differentiating the expression for B(r) and 
equating it with the other expression for B'(r) we obtain 
a single equation relating various derivatives of F, p and 
p. This equation is a non-linear third order differential 
equation for F(r) and due to its length is not explicitly 
shown here. Here we adopt a short-hand notation for the 
equation and write 



mTOV{F, p, p' 



,p",P,p' 



) = Sf, 



(10) 



where S f = Sf(F',F",F"',p,p',p",p,p',p",p"') repre- 
sents a source term. Explicit calculations show that 
in the GR limit, F = 1, the source term vanishes, 
Sf = 0, and Eq. (fTU)) is satisfied, whenever p and p 
satisfy the usual Tolman-Oppcnhcimer-Volkov (TOV) - 
equation. Hence we can view Eq. (TlT)]) as a modified TOV 
equation of metric f(R) theories of modified gravity. 

Given p(r) and p(r), one can solve Eq. (fT0|) for F(r) 
and hence one has a solution for R(r) using the expres- 
sions for B(r) in terms of F. From F(r) and R(r) one can 
determine, at least in principle, F(R) and finally f{R). 
The higher derivative of nature of f(R) theories is appar- 
ent in that Eq. (Tit)]) can have a number of solutions. For 
example, even when matter follows the ordinary TOV- 
equation, one can find non-trivial solutions for F(r) i.e. 
for given standard SSSPF-solution of general relativity 
F(r) — 1 is not the only possible solution. 



B. Boundary conditions 

The outside solution sets the boundary conditions for 
the metric components at the surface of the star, r — 
tq. The field equations are fourth order in A and third 
order in B so that A , A' A%, A'q', B , B' , B%, where 
Aq = A(ro) etc., are fixed. This is in contrast to general 
relativity, where only Aq, A' , Bq are fixed by the outside 
solution. 

In this paper we are interested in solutions which have 
the Schwarzschild-de Sitter -space time (SdS) as the out- 
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side solution. The SdS-metric is 



9iiv = 



( s{r) 



V 





-l/s(r) 









-r 2 

-r 2 



snr 



(11) 



where s(r) = 1 — 2M/r + r 2 /12R and i?o is the corre- 
sponding scalar curvature (in the standard SdS represen- 
tation Rq — — 4A). In the /(i?)-theory framework Rq is 
set by the equation R F(Rq) = 2f(R ). 

This choice of metric fixes the curvature scalar and 
its first derivative at the stellar surface at r : i?(V ) = 
Rq, R'(tq) = 0. The surface of the star is set by the 
conditionsp(ro) — while the energy density is left as a 
free parameter, p(ro) — po- 

Using these boundary conditions in the modified 
Einstein's equations, @, along with the requirement 
RqF(Rq) — 2/(i? ), gives two independent equations: 



Ro 



R"\ = 



s(r) 



dR 2 



R" 



(12) 



Since we are interested in generalized gravity, 
df 2 /dR 2 \ Ro ^ and hence R"(r ) = 0, p = 0. 
Note that this is already a result different from general 
relativity where p can be discontinuous at the surface. 

Derivating the modified Einstein's equations once with 
respect to r and applying the boundary conditions along 
with the new constraints R"(ro) = = po, gives 
R'"{ r o) = 0, p'{r ) = 0. From the continuity equation 
it then straightforwardly follows that also p"(ro) = 0. 
This process can then be continued by derivating the 
modified Einstein's equations once more and substitut- 
ing the boundary values found at previous steps, result- 
ing in three independent equations for the highest order 
derivatives of R, A and B. These can then be solved and 
one can proceed to the next order derivatives. Obviously, 
derivatives of R and A and B are not independent but no 
new information is obtained from their mutual relation. 

Higher derivatives of p and p are not uniquely defined 
and one is free to choose p" and higher derivatives on 
the surface in a way that reflects the desired equation of 
state. 

In summary, requiring that the outside metric is the 
SdS-metric, sets p(r ) = p'{r ) = p(r ) = p'{r ) = 
p"(^o) = on the surface for a general f(R) theory for 
which d 2 f I dR 2 \n a ^ 0. We can already see here how the 
higher-derivative nature of f(R) theories lead to a more 
constrained system than general relativity. 



starts are often approximated by a polytropic equation 
of state, p = up 1 , where k and 7 = 1 + l/n are constants 
and n is often referred to as the polytropic index (see eg. 

The continuity equation, Eq. ([5]), can straightfor- 
wardly be solved for such an equation of state: 



p=«V(l-7)f(A )(l"7)/(27) _ l) 

V B J 



1/(7-1) 



(13) 



where we have used the fact that p vanishes at the sur- 
face. Clearly this condition sets that 7 > 1. Similarly, 
since p' is also vanishing at the boundary, we must re- 
quire that 7 < 2. This is again quite different from GR, 
where 7 is unbounded from above (again see eg. [38[). 



D. SSSPF with a SdS boundary metric 

In light of the discussion presented above, we can now 
consider perfect fluid matter in f(R) gravity theories 
with the SdS metric as a boundary condition. Again, it 
is advantageous to consider a form of the field equations 
where explicit f(R) dependence is eliminated in favor of 
F(r), Eq. ©. Since R'(r ) = R"{r ) = R"'{r Q ) = 0, 
it clear that F'(r), F"(r), F"'(r) vanish on the bound- 
ary ie. F'(r ) = F"(r ) = F'"(r ) = 0. The boundary 
value of F can be solved from the requirement R F(r ) — 
2/(Eo) once f(R) is given, eg. if f(R) = R - p 4 /R, 
Rl = 3p 4 and F = F(Rq) = 4/3. Note that with a 
more general choice of f(R), one can in principle easily 
mimic general relativity ie. Fq — 1 everywhere. How- 
ever, requiring RqFq — 2/(Rq) along with Eq. ©, im- 
plies that p — 3p = everywhere inside the star making 
such constructions physically uninteresting. 

Given a fluid sphere, p(r), p(r) with appropriate 
boundary conditions, one can in principle solve Eq. (|10p 
with the aforementioned boundary values for F(r). Us- 
ing the solution f(R) can then be reconstructed. The 
constant solution, F(r) = Fo = 1, is a solution exactly 
when p(r) and p{r) satisfy the TOV-equation i.e. mat- 
ter is distributed like in general relativity. Vice versa, if 
matter does not obey TOV-equation, the Sf term acts 
as source for the differential equation of F(r) and thus 
the boundary conditions are not strong enough to force 
F be constant, but the solution is more general. An im- 
portant restriction is, however, the boundary conditions 
for p and p making the set of allowed general relativistic 
SSSPF-solutions more restricted in f(R) gravity, eg. the 
standard Shwarzschild fluid sphere with constant density 
is not allowed since po =/= 0. 



C. Polytropic stars 

As a simple example of how the boundary conditions 
already limit the range of allowed solutions, we briefly 
consider polytropic stars. White dwarfs and neutron 



IV. CONCLUSIONS 

In the present letter we have discussed spherically sym- 
metric solutions with non-trivial matter distributions ap- 
plicable to stellar systems in particular. Although we 
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have restricted our analysis to a simplified system de- 
scribed by a perfect fluid, similar conclusions are ex- 
pected to apply in more realistic cases. 

We find that, like in the cosmological case, the dis- 
tribution of matter does not determine the gravitational 
theory uniquely but due to the higher-derivative nature 
of the field equations, different gravitational theories can 
support the same solution. Given the matter distribu- 
tion, p(r), p(r), one can, at least in principle, construct 
a gravitational theory that has the desired solution by 
solving the modified TOV-equation, Eq. (fTU)) . 

By considering configurations that are matched to a 
Schwarzschild-de Sitter metric, we find that such config- 
urations are more tightly constrained than those of gen- 
eral relativity. Again, this is due to the higher derivative 
nature of the metric f(R) theories of gravity that requires 
matching of higher order derivatives at the boundary of 
the fluid sphere than in general relativity. 

As a result, we find that stellar configurations, approx- 
imated by a perfect fluid sphere, can be accommodated 
to an external SdS-solution, whenever p, p and f(R) are 
related by a the modified TOV equation. The conven- 
tional TOV-equations correspond exactly to the choice 



f(R) = R, and departures from the standard TOV- 
equations necessarily lead to a more general gravitational 
action. This phenomenon maybe noteworthy whenever 
modifications to the Einstein-Hilbert action, f(R) = i?, 
are small. Small changes to the EH-action are likely to 
lead only small modifications to stellar models, it. den- 
sity and pressure of the matter may deviate only slightly 
from the ordinary TOV-relation. On this basis it hence 
seems possible that realistic stellar models may be con- 
structed also in f(R) gravity models without violating 
constraints from the solar system. The exact nature and 
whether such solutions correspond to f(R) theories that 
can act as dark energy requires more extensive analysis. 
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APPENDIX 



The set of modified field equations, Eq. ([5]), for a known matter distribution in terms of F(r) reads: 

F 3p 3p Fp 2 Fpp F p 2 FB 1 F' B' F' Fpp' 

" " 272~1 4~ 2r 2 B(p + pf r 2 B(p + pf 2r 2 B(p + pf ~ 2rB 2 2rB~8B 2 r B(p + pf 
Fpp' FB'p' 3F'p' Fp' 2 Fp'p' F" Fpp" Fpp" 



r 



B(p + p) 2 45 2 (p + p) 45 (p + p) B(p + p) 2 2B(p + p) 2 45 2B{p + p) 2 2 B (p + pf 



F F t p t p FB' t F' t ZB'F' Fpp' Fpp' FB'p' 



2r 2 2r 2 B 4 4 2rB 2 2rB 8 B 2 r B (p + pf r B (p + p) 2 45 2 (p + p) 

F'p' Fp' 2 Fp'p' 3F" Fp" 

~45 (p + p) ~ B (p + pf ~ 2B(p + p) 2 ~TB + 2B (p + p) 

F F p p F' B'F' FB'p' F'p' F p' 2 Fp'p' 



2r 2 2r 2 B 4 4 2rB 8 B 2 4B 2 (p + p) 45 {p + p) B(p + pf 2 B (p + pf 
F^__ Fp" 
+ 45 ~ 2B (p + p)' 

This set of equations is not linearly independent and one can solve B(r) and 5'(r)algebraically: 

B'(r) = -2 (2 F 2 p' + r 2 (p + pf (F' - r F") + r F {p 2 + p 2 + p' (F' - r 2 (2 p' + p'j) + p (-(r p') + F" + r 2 p") 
+p(2p-rp' + F" + r 2 p"))) (r 2 (p + p) F' 2 (p + p + r p') + 2 F 2 (p 2 + p 2 - r 2 p' (3p' + p') + r 2 pp" 
+ P (2p + r 2 p"))+rF(-(r 2 F'p'(3p' + p'))+p 2 (3F' - r F") 

+p 2 (3 F' -rF")+rp (-(r p' F") +F'{2p' + r p")) + p (p (6 F' - 2 r F") +r(-{r p' F") + F 1 (2 p' + r p")))j) 
x ((P + Pf (-2 F 2 + r 3 (p + p) F' + r F (r p + r p — F' — r 2 p')) 2 )- 1 
= r rF f2 (p + p + r P ') 2 -2 ^ 6rp' 2 ^ 2rp'p' 2r p" p r (p + p + r pQ F" 

P + P r (p + p) 2 (p + p) 2 P + P P(r)+P 

+F' (-3 + ^4 + rp'(-2p-2p + rp') {p + 2 _ rV ))} 

(p + p) P + P 

x(2F 2 -r 3 (p + p)F' + rF (-(rp) - r p + F' + r 2 p'))^ 1 (1 4) 
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